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We predict that a strong nonreciprocity in the resonance spectra of Dirac quantum dots can be 
induced by the Berry phase. The nonreciprocity arises in relatively weak magnetic fields and is man¬ 
ifest in anomalously large field-induced splittings of quantum dot resonances which are degenerate 
at B = 0 due to time-reversal symmetry. This exotic behavior, which is governed by field-induced 
jumps in the Berry phase of confined electronic states, is unique to quantum dots in Dirac mate¬ 
rials and is absent in conventional quantum dots. The effect is strong for gapless Dirac particles 
and can overwhelm the B-induced orbital and Zeeman splittings. A finite Dirac mass suppresses 
the effect. The nonreciprocity, predicted for generic two-dimensional Dirac materials, is accessible 
through Faraday and Kerr optical rotation measurements and scanning tunneling spectroscopy. 


I. INTRODUCTION 


Quantum dots can be embedded in two-dimensional 
Dirac systems using local gate potentials and point 
charges, as recently demonstrated in graphene [l[ [ 2 [. 
These Dirac quantum dots are defined by nanoscale p- 
n-junction rings, with Klein scattering at the p-n junc¬ 
tions serving as a vehicle for confinement of electronic 
states [sl-Q. Carrier confinement in these ring-shaped 
electron resonators arises due to constructive interference 
of electronic waves scattered at the p-n junction [ 1 , [iq| 
and inward-reflected from the ring. Confined states are 
manifest through resonances appearing periodically in 
scanning tunneling spectroscopy maps[l|, [ 2 [. 

Here we show that this mechanism for electronic con¬ 
finement can be exploited for accessing exotic and po¬ 
tentially useful behavior which is not available in con¬ 
ventional quantum dots. In particular, we predict that 
the Berry phase, a signature topological characteristic of 
Dirac materials Ill4l7j , induces strong nonreciprocity of 
quantum dot resonances in the presence of a weak mag¬ 
netic field B: 


7 ^ ^n,—r, 


( 1 ) 


Here m and n denote the azimuthal and radial quan¬ 
tum numbers, respectively (for optical nonreciprocity, see 
Refs.[l 8 |, Ei). As we will see, resonance splittings of the 
±m states, which are degenerate at H = 0 , grow rapidly 
with magnetic field, approaching values as large as half 
the quantum dot resonance period A^. In particular, 
for the weak B of interest, the effect dominates over the 
H-induced orbital and Zeeman splittings. 


II. SEMICLASSICAL DESCRIPTION 

The Berry phase jumps can be understood from a sim¬ 
ple semiclassical picture describing confined electrons in a 
gapless two-band system. For a confining potential with 
circular symmetry, the resonance spectrum en,m can be 
obtained from the WKB condition for (/^orb = i /c 


the usual orbital phase accumulated along the classical 
path C: 


^orb(^,^) + = 27 r(n + 7), (2) 

where is Berry phase and 7 is a constant [l2[ [3, [ 2 ^ . 
The Dirac band structure, viewed as a Zeeman-type 
Hamiltonian PL = h{p) • tr, where cr = {crx,ay,az) are 
Pauli matrices, gives rise to a geometric gauge field that 
generates the Berry phase, 

PB = ^ dp - {h^\i\/p\h^) = S{C)/2. (3) 

In Eq.(|3]), S{C) denotes the solid angle subtended by the 
vector h = {hx^hy,hz) along a closed path C, and \h±) 
are eigenstates of the two-band Hamiltonian: 


PL\h±)=±\h\\h±). 


(4) 


The Berry phase in a gapless system (hz = 0) can only 
take the values = 0 or Ttt ll, [H, 1 ^. 1 ^. 



FIG. 1. Controlling the Berry phase of confined Dirac elec¬ 
trons using magnetic fields. Shown are semiclassical orbits 
of a massless particle exhibiting topologically distinct orbital 
behavior corresponding to (a) B < Be and (b) B > Be [see 
critical field Be in Eq.([7j)]. The Berry phase, determined by 
the solid angle subtended hy h = {hx,hy,hz) in Eq.(|3]), jumps 
from pb = 0 to pb = 7^ B = Be, see insets [for gapless 
systems hx,y = v (lx,y and = 0, with qx,y the kinetic mo¬ 
mentum (red vectors) and v the Fermi velocity]. Here we used 
m — 1/2, energy e = 1.35/ii’/r*, with r* defined in Ea. fllip . 
BfBe — 0.8 for (a) and BfBe — 1.6 for (b). 
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An external magnetic field can alter the Berry phase of 
the orbits, allowing them to switch between the (/^b = 0 
and ± 7 r types. As illustrated in FigHJ switching can take 
place even in a weak magnetic field. In particular, for 
B = 0 we find (Pb{£, ±m) = 0, whereas for weak nonzero 
fields we find (pB{£,Tn) = ir and = 0 . As a 

result of the tt difference in the WKB condition in Eq. (|2]) 
for the ±m states, the m > 0 and m < 0 families of 
resonances are shifted by half a period, giving rise to a 
large resonance splitting (FigH]): 

£n,m - £n-m - A5/2, (gaplcSs) (5) 

where ^ 10—50meV is the spacing of resonances in 
each family. Equation m describes gapless Dirac band- 
structures, a generalization for gapped systems is dis¬ 
cussed below. 

To illustrate how B controls the Berry phase, we con¬ 
sider a massless particle confined in a radial electrostatic 
potential U{r). This corresponds to h = v{qx,qy,0) in 
Eq.(|4j). In the presence of a uniform magnetic field 5, 
the kinetic momentum q = p — eA is given by 



qr = pr = — U{r)]‘^ /— {mhjr — eBr/ 2)^, 

qg = Pq — eAo = mh/r — eBr {2. 

(6) 

Here v is the electron velocity, and we used the axial 
gauge Ax = —By 1 2^ Ay = Bx/2 to preserve rotational 
symmetry. Because the system is integrable, with con¬ 
stants of motion 5 and m, we can map q to the surface 
of a torus. Figure [2] shows such mapping, with q plotted 
along two curves: Ce in the toroidal direction and Cr in 
the poloidal direction. At a critical B = Be we find that 
the winding number of q along Cr jumps from 0 to 1, 
thus resulting in a 7r-jump of (/9 b • 

The semiclassical quantization of quantum dot reso¬ 
nances can now be obtained from Eq. ([2)) using q in Eq. (| 6 ]) 
evaluated on both C = Cq and C = Cr 0 . This yields 
two quantization conditions for m and 5. For C = Cq, 
Eq.© yields m = -h 70 — (/9B/27r, where epB = tt inde¬ 
pendently of B [see blue curves in panels (b) and (c) of 
Figl2]. Using 70 = 0, we find the anticipated quantiza¬ 
tion of angular momentum m=half-integer. For C = 
instead, we find ^ drpr = 27r(n^ -h 7 r) — where 
ri and r 2 are the classical return points. The half period 
shift in the radial quantization condition results from the 
TT-jump in ifB at H = He. 

While the same semiclassical picture applies to gapped 
Dirac systems {hz 7 ^ 0), there are important differences 
with respect to the gapless case. In particular, the solid 
angle subtended by the vector h{q)^ which now points 
towards the upper hemisphere, is strictly smaller than 
27r; nonreciprocity induced by Berry phase is quenched 
at increasing bandgaps, as will be shown with a more 
detailed quantum model in Figl5l In the limit \hz\ ^ 
\hx,y\, orbital splitting dominates. 


FIG. 2. Topologically distinct mappings of q [Eq.®] to the 
surface of a torus (a), plotted for (b) B < Be and (c) B > 
Be. Indicated with blue(red) arrows is q along the curves 
C 9 {Cr) shown in panel (a), where dotted lines/arrows indicate 
a curve/vector in the bottom surface of the torus. At H = He, 
there is a transition between trivial and non-trivial winding 
of q along Cr. This results in H-induced phase jumps of the 
Berry phase. Here we define g* = s^/v and use the same 
parameter values as in Fig[T] 

The jump in Berry phase corresponds to a transition 
from convex orbits to skipping orbits (Fig(T]). This obser¬ 
vation allows to define the critical field He that induces 
giant nonreciprocity, i.e. the field necessary to reverse 
the electron velocity at the outer classical return point. 
From Eq.® we find qe = rnhjT^iA) — eBeV^iA)I2 = 0, 
with r 2 {s) the outer return point [i.e., qr{r 2 ) = 0]. For 
a quadratic potential model U{r) = /^r^, this condition 
yields 


Bc[T] = 


2hmK 

es 


mK[ey / /im^] 

^[meV] 


(7) 


Using typical values corresponding to recent 
experiments [H, hz ~ 4eV//im^, 5 ~ lOmeV and 

m = 1/2, we find He on the order of 0.3 T. 

Besides the splitting arising at H = He, another 
signature of the nonreciprocal effect is the linear m- 
dependence of He, see Eq. ©• This dependence can be 
understood by noticing that, for larger m, a larger H is 
necessary to induce skipping orbits. As we will see, the m 
dependence of Be gives rise to a peculiar branching pat¬ 
tern of the quantum dot resonances which can be probed 
in spectral measurements away from the quantum dot 
center (see Sec. III). 
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FIG. 3. Magnetic response of quantum dot resonances in a 
gapless Dirac system, (a) The quantum dot is defined by 
the circular p-n ring (dashed lines) induced by a radial elec¬ 
trostatic potential U(r). (b) The magnetic response is dom¬ 
inated by the Berry-phase splitting Asb, which is approxi¬ 
mately half the resonance period Ae. Also indicated in the 
figure is the orbital splitting Asorb- Peak splitting is calcu¬ 
lated from Eq.(|2| for n = 0,1, 2, m = ±1/2, and 7 = 0.6; Be 
is calculated from Eq. 0 ; e, and B^ are defined in Eq. m- 

Importantly, the giant nonreciprocal effect relies on the 
splitting due to Berry phase being dominant over orbital 
and Zeeman splittings. This is the case, in particular, for 
the value Be ^ 0.3 T found in Eq.Q. Indeed, Be is sig¬ 
nificantly lower than the value 5 ll = {As)^/ehv^ ^ IT 
which is necessary for the first Landau level to be larger 
than the resonance period A 5 ~ 25 meV. The strength 
of the nonreciprocal effect is illustrated in FigH] which 
shows the semiclassical spectrum obtained from Eq.® 
for n = 0,1, 2 and m = ±1/2 including both orbital and 
Berry phase splitting. For typical model parameters, the 
splitting A^b ^ Ae/2 induced by the Berry phase jump 
dominates over the conventional orbital splitting A^orb- 
The effect becomes more dramatic at larger n and smaller 
m. Furthermore, the energy ez for the electron Zeeman 
splitting, ez = ^ 10“^ meV, is negligible com¬ 

pared to the characteristic energy of quantum dots (here 
/iB ~ 5.8 • 10“^eV/T is the Bohr magneton). 


with radius smaller than the characteristic length of the 
electrostatic potential, it is legitimate to use a parabolic 
potential model U{r) ^ By using the axial gauge 

Ax = —Byf2, Ay = Bxl2 to preserve rotational symme¬ 
try, the eigenstates of Eq. (jHj) can be expressed using the 
polar decomposition ansatz. 




^imO 


ui{r)e \ 


(9) 


with m a half-integer number. This decomposition allows 
to rewrite Eq.(|8]) as 

/ 5±A/2 dr -\- m/r — Br/2 A A A _ q 

\—dr ^ m/r — Br/2 r ^—5 —A/2 ) \U 2 ) 

( 10 ) 

Here r and B are in units of r* and H*, respectively, 
whereas 5 and A are in units of 5 *, with 

^ 60 nm, 5 * = \/^ 10 meV, 


= {h/e) • {/{K/hvf - 0.2 T. 

( 11 ) 

In these estimates, we considered (gapped) graphene v ^ 
10^ m/s as model system and used a typical value oi n = 
4eV//im^, see estimates below. 

A suitable diagnostics of nonreciprocity, allowing di¬ 
rect access to the quantum dot resonances, is the local 
density of states D{e) inside the quantum dot. Naturally, 
D{e) can be obtained experimentally via the dljdV in 
STS measurements as in Refs. Hi- The quantity B){e^ 
at r = ro can be conveniently written as the sum of m- 
state contributions D{e) = 

Dm{e) = ^{\ua{r = ro)|^)Aa<5(e - £«)• (12) 

a 

Here a labels the radial eigenstates of Eq. (nnii for fixed 
m, and {\ua{r = ro)\'^)\a = 

represents a spatial average of the wavefunction centered 
at r = ro. A gaussian weight is included in the density 
of states to account for the finite size of the tunneling 
region in real STS measurements [l|, [2|. 


III. MICROSCOPIC MODEL 

To supplement the simple semiclassical picture above 
with a microscopic quantum model, we consider the Dirac 
equation describing confined electrons in the presence of 
a uniform magnetic field: 

[va ■ q {A/2)az ± U{r)] 4^(r) = (8) 

Here A is the bandgap and q is the kinematic momentum 
with components Qx.y = —'^^dx^y — eAx^y and qz = 0. This 
corresponds to h = v{qx,qy, A/2v) in Eq.(j4]). Because we 
are interested in eigenstates confined inside the p-n ring. 


Splitting of quantum dot resonances 

Figure m^a) shows the resulting quantum dot spectrum 
as a function of B for gapless Dirac systems, exhibiting 
the H-induced splitting of quantum dot resonances. In 
our calculations, we used ro = 0, Xd/r^ = 0.1, and plot¬ 
ted dD/de in Eq.([T2]) in order to enhance spectral fea¬ 
tures (see Appendix A for details). In agreement with our 
semiclassical interpretation, a half-period splitting is ob¬ 
served in the gapless spectral maps in FigHJ^a). Because 
in FiglU^a) the wavefunction is probed at the center of 
the quantum dot, only small m states {m = ±1/2) con¬ 
tribute to the spectral maps. It is important to stress 
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FIG. 4. Maps of the local density of states as a function of po¬ 
sition ro for a gapless Dirac quantum dot displaying splitting 
of resonances in weak magnetic fields: (a) ro = 0, (b) ro = r*, 
and (c) ro = 2r*. Indicated with dotted lines is Eq.(|7|) for 
half-integer m. The off-centered spectral maps (b)-(c) are 
qualitatively different from the centered case (a) which is sen¬ 
sitive primarily to m = ±1/2 states. Characteristic units for 
magnetic field, is defined in Eq. (HD- Plotted with dot¬ 
ted lines is Eq. © for half-integer m. To enhance spectral 
features, we plot in both panels the derivative of the local 
Density of States in Eq. (HU). 


(a) 



(b) 


Gapped 

< TT 




e/e^ 


e/e* 


EIG. 5. Partial-m contribution to the on-center density of 
states for quantum dots in a) gapless and b) gapped Dirac 
systems. The strong nonreciprocal effect induced by Berry 
phase disappears when a large gap A is opened. As a result, 
resonance splitting is dominated by (a) the Berry phase jump 
in gapless systems, and (b) orbital effects in gapped systems. 
The distinct behavior between (a) and (b) is shown in the 
partial m = 1/2 maps of the density of states [indicated with 
a dotted line is Eq. © with m = 1/2; A/e* = 5 was used in 

(b)]. 


gapped systems is less prominent; in particular, splitting 
is dominated by the orbital contribution. Indeed, the 
peak splitting for the low-energy resonances in gapped 
Dirac systems [e ^ A) can be quantified using a simple 
nonrelativistic model that is valid in the limit A e*. 
In this case, expansion of the Dirac Equation in powers 
of A gives a nonrelativistic Schrodinger Equation for the 
first spinor component 4^i(r): 

[qf^/2A ± t/(r) ± A — e5/2A] 

(13) 

^n,m = ^ (2n ± |m_| ± 1) - 


Here £n,m are the quantized eigenvalues, uj = 
Is. ± e^H^/dA^, and m± = m ±ll2. Here we in¬ 
troduced the orbital magnetic moment = ehv‘^/2 
which can be viewed as an effective Bohr magneton of a 
free Dire particle of mass A/v^. This term is responsible 
for the shift of the resonances with B seen in Figlbfb). 


Self-consistent calculation of the potential profile 


that large m states, which can be probed in off-centered 
STS measurement, are equally susceptible to the Berry 
phase splitting. Figure |ll[b) and (c) show such spec¬ 
tral maps, in which the wavefunctions are probed at (b) 
ro = r* and (c) ro = 2r*. In these cases, there is an over¬ 
lap of peak splitting at different values of H, highlighted 
with fans of Be in Eq.© for varying m (dotted lines). 
At a larger value of ro, states with larger m and 5 can 
be probed. This is indicated by a larger contrast in the 
local density of states induced by such states in FigHJ^b) 
and (c). 

As shown in FiglSl the splitting of the resonances for 


Estimates for n used in Eq.© can be obtained from 
a simple electrostatic model which involves a metallic 
sphere proximal to the graphene plane [Figl6][a)]. This 
model accounts for the fields and charges induced by an 
STM tip on top of graphene, as discussed in Ref.jH. We 
denote with R the metallic sphere radius, and with d 
the sphere-graphene separation. A potential bias SVh be¬ 
tween the sphere and graphene [see Figl6l[b)-(c)] results 
in a spatially varying charge density profile 

Sn(r) «_ + 

47re2(d+r2/2i?)' 


(14) 
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(a) (b) Vacuum 



I SEP'''- 

FIG. 6. (a) Schematics of the electrostatic model, showing a 
metallic sphere of radius R separated a distance d from the 
graphene plane. A potential bias Vb applied on the sphere in¬ 
duces a local variation of the carrier density, different from 
the carrier concentration density Uoo far from the sphere. 
(b,c) Band structure schematics showing band alignment be¬ 
tween the metallic sphere and graphene for (b) large sepa¬ 
ration and (c) close proximity. Here Vcpd is the contact po¬ 
tential difference between graphene and the metallic sphere, 
SVh = Vb — Vcpd, and /io is the Fermi energy under the sphere. 


Here Sn{r) = sgn[/i(r)]/i(r)^/7r(fii;)^ — Uoo is the STM 
tip-induced charge density variation on graphene, with 
the Fermi energy and Uoo the gate-induced carrier 
density far from the tip. Equation ([TT|) is obtained from 
a parallel-plate capacitor model with a slowly varying 
plate separation dc{r) ^ d + r^/2i?. Higher-order terms 
arising due to the curvature of the electric field lines are 
neglected. 

A straight-forward calculation yields a value of k = 
-/i"(0)/2 given by 

K = (15) 

2Rd^/TTW\ 

Here /io is the Fermi energy directly under the sphere, 
and /3 is a dimensionless number defined as 

1 - y/l + 1^1 


For typical values of R ^ 1 /im, d ^ 5nm, SVi, ^ 0.1 V 
and rioo ^ lO^^crn”^, we obtain the value of ^ 4 • 
10“^ eV/nm^. 


(16) 


8e‘^d sgn(/3) 

1 d 

(3 = ^ \e5V\y — Ane^dn 

[hvy ^ 


IV. DISCUSSION 


edge modes which are also present in generic two- 

dimensional structures (2^. 

The large magnitude and tunability of our nonrecip¬ 
rocal effect may help design optical devices, such as 
nanoscale isolators and circulators, which are driven by 
Berry phase. Of special interest are photonic effects 
in Dirac quantum dots. Indeed, electrostatic doping 
can, via the Pauli blocking mechanism, induce a strong 
and tunable electron-photon coupling. This, combined 
with the in situ tunability of the resonance dispersionjj, 
can make Dirac quantum dots useful in miniaturizing 
nanophotonic systems. 

It is instructive to compare our nonreciprocal effect 
with other exotic manifestations of Berry phase predicted 
to occur in Dirac systems, such as Berry phase modifica¬ 
tion to exciton spectra 31, 32|, optical gyrotropy induced 
by Berry’s phase [33| and chiral plasmon in gapped Dirac 
systems [sj, [HI . In realistic electronic systems electron 
decoherence usually hinders observation of such subtle 
effects. We therefore expect that the quantum dot states 
readily available in Dirac materials provide a new and 
optimal setting for locally probing Berry phase physics. 

Since our predictions, such as the strong dependence of 
resonance splitting on A, only rely on the Dirac nature of 
charge carriers, they can be tested in a wide range of ma¬ 
terials. Graphene is the prototypical material to explore 
the case A = 0; graphene on top of axis-aligned hBN 
substrate allows to explore the case A ^ 50 meV [s^, ; 

monolayers of transition metal dichalcogenides such as 
M0S2 allow to explore A of an eV scale 38I-[4 q| . Further¬ 
more, the value of 5* can also be tuned with electrostatic 
potential, as demonstrated in Ref.[l|. 


SUMMARY 

To summarize, quantum dots embedded in Dirac ma¬ 
terials grant access to a new nonreciprocity mechanism 
originating from the Berry phase. This mechanism, 
which is unique to Dirac materials, leads to stronger non¬ 
reciprocity than that for other known mechanisms. The 
anomalous strength of the effect and its in situ tunabil¬ 
ity makes Dirac quantum dots an appealing platform for 
nonreciprocal nanophotonics. The recent introduction of 
Dirac quantum dots in graphene makes these predictions 
easily testable in on-going experiments. 


We point out that the nonreciprocity mechanism con¬ 
sidered above is inherent to Dirac materials. In contrast, 
Faraday and Kerr rotation, two notable examples of non¬ 
reciprocity which can be sizable in two dimensional ma¬ 
terials such as graphene [13, [HI, are also present in gen¬ 
eral semiconductor materials. The same applies to mag- 
netoplasmonic effects, e.g. unidirectional low frequency 
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Appendix A: Computational details 

To solve EQ. (pffi|) above, we use the finite difference 
method in the interval 0 < r < L. The azimuthal quan¬ 
tum numbers are chosen in a finite range, — M < m < M, 
with M large enough to represent accurately the states 
in the energy range of interest. In our calculations, we 
used a system of size L/r* = 10 discretized in N = 600 
lattice sites, with maximum azimuthal quantum number 
M = 31/2. To calculate the density of states, Eq.([T2]), 
we approximate the delta-function 6{e) by a Lorentzian 
6{£) ^ r/7r(5^ + r^). We used a broadening parameter 
= 0.25, and set a Gaussian weight in the spatial 
average . .)x^ of the wavefunction to = 0.1. 
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